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I. INTRODUCTION 


Understanding the environment has been a problem mankind has dealt with for 
ages. Our knowledge of the earth was initially limited to that which we could see 
with our own eyes. Gradually our ability to sense our environment grew as new 
instruments were developed. When space satellites came onto the scene, a rich 
resource of information concerning the world around us suddenly opened. From visual 
and infrared images of the earth, the truly dynamic nature of the planet we live on 
could finally be appreciated. 

The years following the first attempts at retrieving environmental data from space 
have shown advances that parallel those made in computer technology. Devices have 
been implemented which have made use of virtually all of the electromagnetic (EM) 
spectrum. Each remote sensing device was designed for the purpose of detecting some 
form of radiation either emitted or reflected by the earth. To observe particular 
features of the atmosphere or ocean, methods were developed which would utilize 
particular frequencies to penetrate through obstacles in the path of this radiation and 
arrive at its target parameter. 

Wind speed over the ocean has proven to be a particularly difficult parameter 
to measure, but most worthwhile to pursue. Vast areas of ocean have been ‘data 
deserts" for years. Remote sensing of these regions using microwave frequencies can 
supply data to improve ‘data starved" numerical input for environmental models. 
This, in tum, should offer a better understanding of synoptic weather systems, 


improving on the visible and infrared technologies currently used today. Disseminating 


information concerning the global surface wind field to all vessels, would enhance 
preparations for oceanic transit. Ships already underway would be better prepared for 
heavy winds and seas. 

In the case of wind speed over the ocean, two separate techniques have been 
developed. The first is the active radar device known as the GEOSAT altimeter which 
utilizes pulses of EM energy at 13.7 GHz to analyze the surface of the ocean. Though 
its prime objective was to map the earth’s geoid, its signal can also be used to measure 
surface wave height and wind speed over the ocean. The second device known as the 
SPECIAL SENSOR MICROWAVE IMAGER (SSM/I) utilizes a totally different 
approach to the problem of wind speed measurement. Instead of transmitting its own 
radiation to the earth and then receiving a returned signal for analysis, it passively 
measures upwelling radiation from the earth’s atmosphere and oceans. 

How these two methodologies compare, is the subject of this thesis. The data 
utilized in this comparison represents several areas of the world where both satellites 
were co-located within one hour of each other’s orbit. In the pages to follow, details 
of how each system operates and the physics behind the analysis of the data, will be 
explored. Correlation studies made for each region will disclose similarities and 
differences in the wind speed measurements. A particularly interesting case will be 
presented where several environmental parameters have interfered with the sensors 
ability to accurately access the surface wind speed. 

Chapter two details the basic physics applicable to the SSM/I and GEOSAT 
remote sensors, as Well as the operation of both systems. Chapter three will discuss 


ground truth verification studies made for each system. Here, individual strengths and 


weaknesses are taken into account prior to the comparison study. The experiment is 


presented in chapter four, followed by conclusions and recommendations in chapter 


five. 


we 


Il. FUNDAMENTALS OF MICROWAVE REMOTE SENSING 


The basic physics behind the remote sensing techniques examined in this thesis, 
will be discussed in this chapter. The first section will cover the propagation 
characteristics of EM radiation. This will be followed by an overview of the index 
of refraction, blackbody radiation, emissivity, and the radiative-transfer equation. This 
should provide an ample background for understanding the operation of both the SSM/i 


and the GEOSAT altimeter. 


A. PROPERTIES OF EM WAVE PROPAGATION AND ATTENUATION 


Microwave remote sensing depends on the modification of EM radiation by the 
environment. Changes take place in the phase and intensity of an EM wave as it 
passes through the atmosphere and interacts with the ocean surface. These alterations 
of the wave’s character are interpreted by the remote sensing device as variations in 
the state of the environment. We will first describe the EM wave equation and then 
develop the relationships that govern the changes that take place in the wave as it 
passes through matter. 

EM radiation consists of alternating electric and magnetic fields [Ref. 1], which 
are classified by their frequency of fluctuation. A single EM wave vector, E(x.y,z), 
can be broken down into its individual components, E,, E,, and E, for a simplified 


analysis (Fig.l). The propagation characteristics of a single component of an EM wave 






Direction of 
Propagation 


Figure I Electromagnetic wave vector and components. 


can be described by the following equations: 


E, = E, chin) 
k* = @’pe - imp 
@ = 2nf 
B= HoH, 


E = €:60Scexe - 1€”) 


tT] 
lM 


component of the electric field intensity 


ee 
IM 


initial electric field magnitude 
k = complex wave number 
@ = frequency in radians per second 


f = frequency in hertz 


(1) 
(2) 
(3) 
(4) 
(5) 


lt, = permeability of free space (4m x 10’Hz/m) 


relative permeability of a matenal (=1) 


H, 
€ = permittivity of a medium 

€, = permittivity of free space (8.85 x 10"? Farads/m) 

e’.e” = real and imarginary components of the dielectric constant 


For most cases in remote sensing applications, the value of pp, = 1, so that 


equation (2) can be written as: 


kK? = w'p.ele’ - ioKe,0)] (6) 
where: 
€ = [e’ - io/(€,@)] the complex dielectric constant 
oO = conductivity of the medium 


In free space, i.e., a vacuum, €, = 1, o = 0, and the EM wave velocity is given 
as: 
c. = 1KM.€,) (7) 
For EM waves traveling through matter, such as the ocean or atmosphere, the velocity 
is determined through the equation: 
a= ak =e (8) 
From (7) and (8), we can define the index of refraction n, which is related to €, by 
the formula: 
n=c/c = Ve. (9) 
Equation (9) is umportant for understanding how the EM wave is effected as it passes 
through the atmosphere and interacts with the ocean below. Since n is related to the 


complex dielectric constant €,, radiation passing through matter will experience a 


modification in speed and intensity. The wave will be slowed due to the resistance 
encountered by an induced electric field produced in this medium and attenuated 
through the mmaginary component of n. To understand how the imaginary component 
can attenuate the wave, n must be examined in more detail. The index of refraction 
written in its component form: 

n=a+ if (10) 


where & and 6 represent the real and imaginary parts of n and are written as: 


o? = (le,| + ©’/2 le) 
Bv= (leah = © 2 (12) 
le, || = fe? + o/(e,0)"] (13) 


Returning to equation (8) for the wave number k and substituting the component form 
of n into that expression, we can rewrite k as: 

k = wma/c + iwB/c (14) 
By replacing the expression for k in equation (14) into the original expression for the 
propagating wave, equation (1), we have: 


E = E. el-itoaBre - O)] eloprie ] (15) 


Me 
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where = -@/c is known as the attenuation coefficient. With a negative value in the 
exponential, this component acts to reduce the value of the electric field as the EM 
wave passes through matter. The distance the radiation travels through some medium, 
in which the intensity of the wave is reduced by a factor of exp(-1), 1s known as the 
skin depth d and 1s given by : 

d= "ce (@ummor "d = c/( 27115) (16) 


For a frequency of 5 GHz, the value of d is 0.5 centimeters. The SSM/I operates 
between 19 and 37 GHz and the GEOSAT altimeter operates at 13.5 GHz, and so the 


skin depths for the instruments in this study are on the order of millimeters. 


B. POLARIZATION 


In addition to the variation of the intensity and frequency of the wave, microwave 
radiation is also characterized by its polarization. Polarization is defined as the 
orientation of the electric field vector relative to the surface it is reflected or 
emitted [Ref. 2]. Polarized radiation is classified as either horizontal, in which case 
the field is parallel to the surface, or vertical, whereby the field is perpendicular with 
respect to the surface. A change in surface conditions and viewing angle, wil affect 
the vertically and horizontally polarized EM waves in different ways. The amount of 
radiation emitted from a surface (emissivity) depends on how the radiation is polarized 
(Fig.2). The SSM/I’s sensors can detect both horizontally and vertically polarized 
radiation. Certain combinations of polarized frequencies are used to allow particular 
environmental parameters to be measured separately with a minimal of interference 
from the other parameters. 

Assuming a smooth (specular) surface, the reflection of the spectral radiance, 
R, from such a surface is written in equation form as: 

R(A;8,,0,-7) = L(A;8,.,)/L,(A;6,.0,-10) (17) 
where: 


r 


wavelength of the radiation 


8,0 = spherical coordinate angles 
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Figure 2 Polarized microwave radiation at 3 cm wavelength as a 
function of emissivity and incidence angle. 


LL, = reflected and incident radiance respectively 


The spectral reflectance of a plane dielectric surface is determined through the Fresnel 
reflection coefficients as follows [Ref. 3): 
R,, = [(p - cos8)? + q’]/[(p + cos6)? + q? (18) 
Ry = [(€'cos® - p)? + (€”cosO + q)’}/[(e’cos8+ p)? + (€”cos®@ + q)’] (19) 
where R, and Ry, are the vertically and horizontally polarized components of the EM 
wave, and € = €’ - i€”; the complex dielectric constant of the surface. The quantities 
p and q are: 


p = (1/V2)([(e’ - sin’)? +e”7]" + [e’ - sin?0])” (20) 


(1/N2){[(e’ - sin?@)? +e)" + [e’ - sin?@]}* (20) 


p 


q = (IN2){{(e’ - sin’@)? +e”7J" - [e’ - sin?@])" (21) 
where @ is the incidence angle. 
In the case when 9 = Q, such as with the viewing angle of the GEOSAT 
altimeter, R,, and Ry are reduced to the form: 
Ry = Ry = [(p - 1) + qQ’[(p + 1)° + gq’) (22) 
= [(a - 1)’ + BY/{(a + 1) + BY 
where the complex index of refraction is given by equation (10). Since n’ = e’ - ie”, 
the Fresnel reflection coefficient at normal incidence — can take the form: 
Ry = I(n- Dm + 1)/? (23) 
The spectral emissivity of a surface is related to the spectral reflectivity through 
Kirchoff’s law (20). We can relate the spectral reflectance R to the emissivity in 
both horizontally and vertically polarized forms: 
e,(8) = 1 - R,(8) (24) 


e,(8) = 1 - R,(0) (25) 


C. BRIGHTNESS TEMPERATURE 


The purpose of this section is to develop the relationship between the radiation 
emitted by the ocean’s surface and the upwelling temperature detected by the SSM/I 
known as the brightness temperature, T,. All objects, except those existing at absolute 
zero, emit radiation. An object is classified as a perfect emitter, known as a black 
body, if it radiates energy at the maximum rate predicted by Planck’s radiation law. 


As the temperature of an object increases, the total amount of energy emitted also 
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increases. Most objects, however, emit at a rate less than that of a blackbody. The 
term “grey body’, is used to describe an object which emits radiation at a reduced rate. 
Variations in the temperature of a grey body relative to the assumed blackbody 
temperature of the earth can be identified and attributed to changes in the environment, 
such as wind speed, water vapor, and precipitation. 

Terms commonly used in discussing brightness temperature, or radiation 
measurement in general, are radiant energy, radiant flux, irradiance, and emittance. 
The term radiant energy, Q, refers to the amount of energy carried by the radiation 
(Joules). The radiant flux, ® = dQ/dt, is a measure of rate of energy transport 
(Joules/sec). This flux defines the power available to produce a response in a 
measuring device. The irradiance, E = d®/dm’, is used in reference to the radiance 
incident on a surface. The radiance emitted from a surface is known as the emittance, 
M = d@®/dt. The character of the material from which radiation is either emitted or 
reflected is known as its emissivity. Each surface has a particular value, for example, 
the emissivity of salt water is 0.3, while that of dry soil is 0.92. It 1s this difference 
in emissivity that allows the SSM/I to distinguish between land and sea surfaces. 

The GEOSAT relies on the degree of reflectivity to detect changes in the ocean's 
surface. The GEOSAT and SSM/I observations are linked by Kirchoff’s law, which 
relates emissivity and reflectivity in the following manner: 

e=l-r (26) 
where r is the reflectivity of the surface and e is the emissivity. 

The ratio of spectral brightness emitted in a direction (9,0) by a real body 


relative to that which is radiated by a blackbody is the spectral directional emissivity 


1} 


e of that body, where € is a function of wavelength A and spherical angles 6 and 8. 
The microwave emission from the ocean surface is given by eB,, where B, follows 
the Rayleigh-Jeans approximation of Planck’s Law. Planck’s Law states that the spectral 
brightness per unit of frequency is expressed as: 

B, = @hv (epee [eo ay) (27) 
where: 


k = 1.38 x 10” J/K, Boltzman’s constant 


h = 6.63 x 10% J s, Planck’s constant 
T = temperature in degrees Kelvin 
Vv = frequency of radiation in hertz 


c = 3 x 10° m/s, speed of light in a vacuum 
The Rayleigh-Jeans approximation is made at microwave frequencies, where hv/(kT) 
<< 1, so that: 

een s= | eeeiyy (ke) 
Substituting the above relationship into equation [18] gives 
Be (2y Keon (28) 

Equation [19] states that the spectral brightness is proportional to the temperature of 
the body emitting radiation. The radiance received by a microwave sensor is known 
as the brightness temperature. The brightness temperature T,, is given by €T and will 
be discussed in a later chapter. 

The following section will cover the topics of absorption and the interaction of 
radiation with surfaces. These factors contribute to the variation in brightness 


temperatures as measured from satellites. 
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D ABSORPTION AND SCATTERING OF RADIATION 


Radiation emitted by the ocean surface must pass through the atmosphere once 
before being detected by the SSM/I, and twice for the GEOSAT altimeter. 
Consequently, interactions which take place between the atmosphere and the radiation 
must be understood in order to interpret a retrieved signal. In the atmosphere, energy 
is both scattered and absorbed. For microwave frequencies, the primary absorber is 
water vapor, while scattering is predominantly due to rain. Any instrument used to 
measure wind speed over the ocean must utilize frequencies where the effects of 
atmospheric absorption and scattering are small. The GEOSAT altimeter frequency was 
chosen such that atmospheric effects would be small. In addition to observing the 
ocean surface, the SSM/I measures atmospheric parameters such as water vapor and 
precipitation. This allows the SSM/I to adjust for atmospheric parameters which would 
interfere with the measurement of wind speed. Questionable wind speeds, based on the 
presence of these interfering parameters, are identified at ground stations when the data 


is analyzed. 


KE. THE RELATION OF SURFACE ROUGHNESS TO EMITTED RADIATION 


The roughness of the ocean causes the reflected and emitted EM waves to be 
sent out in numerous directions. Remote sensors of microwave radiation must be able 
to account for this variation in directionality, hence the viewing angle becomes an 


important factor. The received radiation 1s also a function of the dielectric constant 


13 


of the surface. These properties can be related if we, (1) assume that the surface is 
made up of many small facets, (2) examine only a single facet, or plane dielectric 
surface and apply the Fresnel Reflection Coefficients to that surface to obtain the 
spectral reflectance of that surface. Finally, (3) integrate over all the facets in the area 
observed by the sensors to obtain the spectral reflectance for that surface. 

The solution to the Fresnel equations requires that either the complex dielectric 
constant, €,, or the index of refraction of the surface as a function of frequency be 
known. In the case of the ocean, the dielectric constant is a well known parameter 
through the work of Lane and Saxon [Ref. 4]. The dielectric constant is represented 
by the Debye equation: 


€ =e, + ((€€.)/1 + (iwt)’* - 16/we,] (29) 


where: 
@® = 2nf, radiation frequency 
€., = dielectric constant at infinite frequency 


the static dielectric constant 


mM 
II 


t = the relaxation time 


the ionic conductivity 


Oo 


Q = an empirical constant 
The parameters €,, t, and © are all functions of the material which is reflecting or 
emitting the radiation to be received by microwave sensors. The conductivity 6 of 


the sea water increases the complex portion of the dielectric constant, however, it is 


important to note that the influence of sea salts on the emissivity of the ocean is small 
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for frequencies greater than 5 GHZ. Since both the SSM/I and GEOSAT sensors 
operate above 5 GHZ, the salinity of the sea water has no effect. 

The above description for spectral reflectance and emittance from a single facet 
can now be applied to the entire rough surface. This is accomplished by integrating 


over all the facets and taking a statistical average. 


F. THE RADIATIVE TRANSFER EQUATION 


The radiative transfer equation contains terms which relate radiant energy from 
the surface and atmosphere to the radiant energy observed by a satellite. Before 
energy can reach a Satellite sensor from the ocean surface, the upwelling radiation is 
attenuated and joined by additional emission from the atmosphere, so that the observed 
temperature detected by the SSM/I is given by: 

T = T,(1 - t) + tleT, +(1 - e)(1 - t)T,] (30a) 
or 


T = et[T, - T,(1-t)} + T,( - 2) (30b) 


where, T, is the effective temperature of the atmosphere, t = e‘’, Tt 1s the optical 
depth, T, is the sea surface temperature, and e is the emissivity. The brightness 
temperature of the radiation incident on the ocean surface by the atmosphere is 

T,(1-t). The amount that is reflected by the sea surface is written as (1 - e)T,(1 - 
t), while the amount emitted is eT,. It follows that the total upwelling brightness 
temperature from the sea surface is eT, + (1 - €)(1 - t)T,. The signal from the 


surface is €T,. The other terms are considered noise with regard to surface wind 
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speed analysis, but are useful when detecting water vapor, liquid water, and 
precipitation. If the optical depth t = 0, then the atmosphere is transparent and allows 
just the brightness temperature of the ocean surface to be seen by the sensor. Equation 
(30b) would then be reduced to simply T = eT,. This the ideal case for measuring wind 
speed over the ocean. If, on the other hand, t >> OQ, then T = T,, te., the SSM/I 
would be looking at the environmental changes in the atmosphere only. The 
frequencies and the method used by the SSMI and GEOSAT to minimize atmospheric 


interference will be discussed in sections devoted to instrument operation. 


G. THE TRANSFER OF MOMENTUM FROM THE ATMOSPHERIC SURFACE 


LAYER TO THE SEA SURFACE 


Waves and foam generated by the wind contribute to the roughness of the sea 
surface and effect both microwave emission as measured by the SSM/I and microwave 
scattering as measured by the GEOSAT altimeter. The roughness has the effect of 
altering the incidence angle of radiation being emitted or reflected and changing the 
emissivity of the sea surface itself. The roughness of the sea surface is an indication 
of the transfer of momentum into the sea from the turbulent atmospheric surface layer. 
An extrapolation process 1s performed to relate the wind speed in the atmospheric 
surface layer to the sea surface roughness. This section will briefly describe the 
transfer of momentum from the atmospheric surface layer to the sea surface. 

Businger [Ref. 5] describes the formulation of the turbulent transfer of momentum 


from the atmospheric surface layer to the sea surface as a semi-empirical relationship. 


16 


The appropriate equation can be developed by utilizing the turbulent transfer coefficient 
K in the equation for momentum transfer, where: 

F, = -pKdU/ez (31) 
Here, p is the density and U is the mean horizontal wind speed. Businger defines 
the atmospheric surface layer as the lower 20 meters of the atmosphere. Here the 
interaction with the surface is strong and the layer adjusts to the surface conditions 
rapidly. We can therefore assume a quasi-steady-state. The momentum transferred 
through to the ocean surface exerts a drag force per unit area on that surface known 


as the surface stress t. By Newton’s third law, the momentum equation is given as: 





F, =puw = -t (32) 
The friction velocity u. is commonly used as a scaling velocity, where: 
u. = (t/p)” (33) 
Equation (31) can now be written as: 
K_dU/dz = u.? (34) 
where u. is independent of height. To generate a wind profile from the surface to a 
height of 20 meters, an assumption must be made about K,. By dimensional analysis, 
K,, has the dimensions of velocity times length. We use u. to scale the velocity and 
z to scale the length. Therefore, we assume that: 
Ker= Kinz (5) 
where, k is a constant of proportionality, the von Karmann constant. Combining 
equations (34) and (35), we have: 


dU/dz = u./(kz) (36) 


i, 


To avoid an infinite shear at the surface, when z = Q, a finite surface roughness z, is 
introduced into the equation. This now becomes: 

dU/dz = u(k{z + z.]) (37) 
Upon integration we have the form: 

Uju. = 1/k Inf(z + z,)/2.] (38) 
which is the logarithmic profile for the wind from the sea surface through the 
atmospheric surface in a neutrally stable atmosphere. This is the common profile used 
to extrapolate surface wind speed to a height of 10 or 20 meters, which are the typical 
heights at which sea buoys and ships measure wind speed at sea. The neutral case is 
used primarily for its simplicity and also for its adequate approximation of the surface 
layer wind field. Figure (3) [Ref. 6] represents a semi-log plot of wind speed from 
the surface through the atmospheric boundary layer for neutral, stable and unstable 
conditions. It would be impractical to measure the stability of the atmosphere for all 
wind speed measurements. Fortunately, the neutral profile of wind speed approximates 
the atmosphere in most situations. 


Fig. 3 is derived from setting equation (37) equal to a wind shear function, ®,,, 


such that: ie Ge Se Re 
where: 
ol te (ea ez: for z/L > O (stable) (39a) 
(Dei for z/L = O (neutral) (39b) 
De = |i-+ [to7 bee for z/L < O (unstable) (39c) 


and L is the Obukhov length. 
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Figure 3 Typical wind speed profiles versus static 
stability in the surface layer. 
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Ill. MEASURING WIND SPEED OVER THE OCEAN WITH THE SSM1 


The SSM/I is a microwave radiometer which measures upwelling radiation from 
the earth’s surface at five frequencies and two polarizations. Using frequencies 
associated with atmospheric windows, the SSM/I measures upwelling thermal 
microwave radiation from the ocean surface. This radiation is called the brightness 
temperature. From equation (30b), where T, = eT, for an atmosphere with 100% 
transmittance, the brightness temperature is a function of the emissivity, e, and 
temperature of the sea surface, T,. The sea surface temperature is relatively constant 
compared to the change in emissivity. An increase in T, is thus caused by an increase 
in the sea surface emissivity. A change in surface roughness and thus wind speed, wil 
change the emissivity of the surface. Thus, a particular wind speed at the sea surface 
can be related to a certain brightness temperature. 

The SSM/I utilizes the 19.35H, 22.235V, 37H, and 37V GHz channels for 
measuring wind speed, where H and V represent horizontal and vertical polarizations 
respectively. The 19.35H, 37H, and 37V GHz channels are "windows" in which there 
is nearly 100 % transmission of upwelling radiation through the atmosphere. The 
22.235 GHz channel is in the center of a water vapor absorption line and is used to 
measure water vapor concentration. Detection of water vapor is vital to the 
determination of the validity of the wind speed detected, since the 19.35H, 37V, and 
37H channels are subject to some small amount of atmospheric interference which 


require that an atmospheric correction be applied. 
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The algorithm used by the SSM/I for measuring wind speed assumes the 
1uughness of the ocean surface can be treated as a conglomeration of many tiny 
emitting planes whose angle of inclination is altered by the force of the wind. The 
variation in the sea surface slope is used in calculating the emissivity of the surface 
by assuming a Gaussian slope distribution of the angles of the slope. A change in 
the distribution of these angles corresponds to a change in the root mean square slope 
of the sea surface and distribution of the emissivities associated with the variance of 
sea surface slopes. The resulting change in emissivity causes a change in brightness 
temperature through the Fresnel equations (Eqns.22 and 23). The equation used to 
describe the dependence of derived wind speed on the surface slope variance, 0° 


[Ref. 7] , is as follows: 


o°(f) = (0.3 + 0.02f) (0.003 + 0.48U,,) for f < 35 GHz (40a) 


o’(f) = 0.003 + 0.48U,, for f = 35 GHz (40b) 


Here, f is defined as the frequency in GHz and U,, as the wind speed measured in 
meters/second at a 20 meter height. 

Foam generated by breaking waves occurs at wind speeds greater than 7 m/s. 
The affect of foam on the microwave emission from the sea surface is not well 
understood and there are varying theories in literature. Swift [Ref. 8] cites an 
experiment over the Salton sea using a 19 GHz radiometer to measure wind speed 
over a foam covered sea. Swift states that the results were tens of degrees higher than 


theory and could be attributed to the presence of foam. Wilheit [Ref. 9] treated foam 
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as a relatively non-reflective surface and eliminated it from the radiative transfer 
equation. The reflectivity that does occur comes from the dielectric interface around 
the foam bubbles. Hollinger et al [Ref. 10] treat foam as “a substance which partially 
obscures the ocean surface in a manner independent of polarization but dependent upon 
frequency’. He defines the foam fraction, F,, as 
F, = 0.006 (U,, - 7) (1 - e”° ) for U,, > 7 m/sec (41a) 
F,= 0 for U, < 7 m/sec (41b) 
This equation is developed from the observation that the sea surface brightness 
temperature increases linearly for wind speeds in excess of 7 m/sec, and that there is 
an absence of foam development for speeds less than or equal to 7 m/s. The foam 
fraction is also used in the empirically derived SSM/I wind speed algorithm to account 


for the change in emissivity that occurs with the presence of foam. 
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IV. MEASURING WIND SPEED OVER THE OCEAN WITH THE GEOSAT ALTIMETER 


The primary mission of the GEOSAT satellite was not wind measurement, but 
measurement of the geoid and sea surface height. However, wind is computed as a 
by-product and with an accuracy of +/- 2 m/s [Ref. 10]. To measure wind speed, the 
GEOSAT satellite utilizes an active radar operating at a frequency of 13.9 GHz called 
an altimeter. This instrument transmits a pulse of microwave energy to the ocean 
surface and measures the return signal. The altimeter [Ref. 11] uses an automatic gain 
control (AGC) loop to normalize the amplitude of the return signal from the sea 
surface. The AGC is a measure of the backscatter cross section of the sea surface, 0 , 
and is related to the wind speed. The greater the wind speed, the lower the cross 
section, thus the weaker the return signal from the ocean surface. The backscatter 


cross section 1s obtained from the following equations [Ref. 12]: 


AGC [dBm] = 10 In (P//P,) (43) 
where: 
P, = the received power 
P, = the reference power (1 mW) 


The received power P,, is defined in terms of antenna gain G, transmitted power P,, 
the operating wavelength A, the distance to the sea surface R, and the scattering cross 
section 0 . The relation appears as follows: 


P= [A*(4n)] P,Go/ R‘dA (44) 


Ze 


Assuming that the sea surface slope is both gaussian and isotropic in its probability 
distribution, we can define 0 as follows: 
o =[ IR(O)|? / S? ] sec’ exp [ -tan’O / S? ] (45) 
where: 
© = angle of incidence 


R(O) = Fresnel reflection coefficient for the air-sea interface at normal incidence 


The Cox and Munk [Ref. 13] relationship for average upwind/crosswind mean 
Square sea surface slope 0° and wind speed is given as: 
o=a+BU (46) 
where, U is the wind speed in m/s and a, B are empirical constants of proportionality 
which depend on the probability density of the sea surface slope. By making the 
assumption of an isotropic Gaussian surface, the separate effects of the upwind 
component and the crosswind component are ignored [Ref. 14]. By ignoring wind 
direction when deriving a relationship between wind speed and surface slope, error in 
the surface wind speed will be introduced. 
If we let © = 0, since we are observing at the nadir angle, the equation for 6 
becomes: 
o = IR(O)I7*(a + BU) (47) 


The equation for the received power is 


P. = (A*/(4n)’) P, (G*/h*) o A (48) 
where: 
h = Satellite altitude 
A 


This relationship is used in the Brown algorithm for the actual calculation of wind 


area illuminated by the altimeter 
speed [Ref. 15] at the Fleet Numerical Oceanographic Center. Brown’s algorithm 1s 


a result of a comparison of wind speed data received from sea buoys to normalized 


backscattered cross sections derived from altimeter power measured on board the 
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GEOS-3 satellite and applied to GEOSAT. The algorithm is presented in two stages. 
The first section 1s as follows: 


W, — exp Kao x o°/10) _ B)/A] (49) 


where: 


A = 0.080074 B = -0.124651 for o < 10.12 dB 
A = 0.039893 B = -0.031996 for 10.12 < o < 10.90 dB 
A = 0.015950 B= 0.017215 for o = 10.9 dB 


In the above equation, the backscattered cross section is given as decibels relative to 


1 m’. 
A second equation 1s used to complete the derivation of surface wind speed. This 
was necessary to adjust for the skewness in the distribution of the differences between 


the buoy and altimeter wind speed measurements. The second part of the algorithm is 


as follows: 
Ui. = 2 a, W, , for W, < 16 meters per second (50a) 
Ce ee for W, > 16 meters per second (S0b) 
where: 
a, = 2.08779900000 
a, = -0.36499280000 
a, = 0.04062421000 
a, = -0.00190495200 
a, = 0.00003288189 


U,, 1s the wind speed at a height of 10 meters. This algorithm appears to measure 


the wind to within +/- 2 m/s. 
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V. GROUND TRUTH STUDIES 


A. GEOSAT Wind Speed Ground Truth Study 
Shuhy [Ref. 16] made a ground truth study for the wind speeds obtained by 


GEOSAT using National Data Buoy Center buoy data. To avoid possible interference 
from land, most of the buoys used were located in deep water away from any 
coastlines. Shuhy chose to restrict the comparison to satellite passes within 80 
kilometers and 90 minutes of a ground truth point. The data was classified as either 


excellent, good, or fair based on the following criteria: 


DISTANCE FROM TIME DIFFERENCE 
CLASSIFICATION GROUND TRUTH IN OBSERVATIONS 
excellent 25 KM 30 min 
good 25 - 60 KM 60 min 
fair 60 - 80 KM 90 min 


No data outside 80 kilometers or 90 minutes from crossing were used. The results of 
the comparison are seen in Fig. 4. Shuhy concluded from regression analysis that the 
sub-categorization described above yielded no significant differences. A total of 256 
data points were used with an mms error of 1.54 m/s. The comparison demonstrated that 
GEOSAT wind speeds are higher than the ground truth for speeds less than 6 m/s and 
for speeds greater than 6 m/s, GEOSAT winds are lower than those obtained on the 


surface. Shuhy notes, however, that there is significant scatter in the data and that 
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Figure 4 Geosat vs ground truth winds at a height of 10 meters. 


spatial variability must be examined to determine which comparisons might be invalid 
due to high variability between buoy and satellite wind speeds along the satellite track. 

When ground truth data is compared to off-nadir GEOSAT data, there appears to 
be a correlation between the magnitude of the angle and how well the GEOSAT and 
buoy wind speed compared. For angles greater than 1.1 degrees, the correlation was 
poor. Angles less than 0.4 degrees off nadir exhibited good correlation in the 


GEOSAT/buoy wind speed comparisons. In this study, the nadir angle was not 


pe 


available as the GEOSAT data used was preprocessed by Naval Ocean Research and 
Development Activity (NORDA), however all data with off nadir angles greater than 


1° is not transmitted to FNOC. 
B. SSM/I WIND SPEED GROUND TRUTH STUDY 


The SSM/I algorithms have undergone four revisions since the system was put 
into space in late 1987. The study completed by Hollinger and Poe (Ref. 18] includes 
600 SSM/I and buoy data sets. Comparisons of wind speed measurements were made 


within 30 minutes and 25 kilometers of each other (Fig. 5). 
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Figure 5 SSMI wind speed comparison with ground truth from moored 
buoys (Hollinger and Poe). 
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The SSM/I winds are extracted from a data set by the "D-matrix" method. This 
is an empirical/statistical approach that chooses the most probable atmospheric and 
surface conditions that would produce a characteristic set of measured brightness 
temperatures. The "D-matrix" contains the reference brightness temperatures for the 
SSM/I. The SSM/I retrieved brightness temperatures are compared to the "D-matrix" 
and from this comparison the ocean and atmospheric conditions that were observed 
by the SSM/I are determined. The overall scatter of data points represents a +\- 2 
m/s standard deviation between 3 and 15 m/s. At this time, there is insufficient data 
for wind speeds greater than 15 m/s to make an accurate comparison of SSM/I and 
buoy wind speeds. 

Interfering factors that could reduce the accuracy of the analysis are heavy 
concentrations of water vapor and the presence of cloud water and rain. In this study, 
areas exhibiting these conditions were included in the data to observe their influence 
on retrieved wind speeds. 

Positioning error was discovered by as much as 20 to 30 km for shift in pixel 
latitude and longitude coordinates. However, Hollinger and Poe state that this did not 


interfere with the overall] validation of the wind speed algorithm. 
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VI. SSM/T AND GEOSAT ALTIMETER WIND SPEED COMPARISON 


A. DATA ACQUISITION 

The data base used was acquired from Fleet Numerical Oceanography Center 
(FNOC) in Monterey, California. All data obtained for this study was gathered in 
near-real time using a personal computer (PC) with a 1200 baud modem. Co-located 
passes of the DSMP and GEOSAT satellites within one hour of each other were 
required. The position and time of ascending node for the orbits of each satellite 
were calculated well in advance. These orbit tables were updated every two weeks to 
account for errors in calculating the exact times of ascending node crossings. 
Adjustments had to be made to the GEOSAT orbit table to account for the periodic 
maneuvering of the satellite that would keep it in a 17 day repeat orbit pattern. To 
avoid any microwave interference from land stations, areas over the Open sea were 
chosen. 

Several problems were encountered in obtaining data in this manner. The average 
length of time to download a single SSM/I data set was at 1200 baud was 
approximately 0.5 to 1 hour, depending on its size. Near real time SSM/I and 
GEOSAT data at FNOC is accessible for 6 hours or less before it is removed from 
memory. As a result. many data sets obtained were incomplete as the data had been 
removed during the downloading process. I was unable to determine at what times a 
particular orbit would be resident at FNOC and for what length of time it would be 


resident. I consequently worked with the data available at the trme. Another difficulty 
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arose in the form of the GEOSAT data format. The GEOSAT data obtained from 
FNOC was pre-processed by NORDA in Bay St. Louis, Mississippi. The format of 
the GEOSAT data base included only the latitude, longitude, wind speed, and 
significant wave height. Raw data would have been ideal, so that various wind speed 
algorithms could have been tested. Occasionally incomplete data sets were obtained 
where wind speeds were missing for large areas of an orbit. This was due to data 
transmission problems occurring between NORDA and FNOC. Despite these 
difficulties, five data sets were obtained from three different regions of the globe: the 


North Atlantic, the Gulf of Alaska, and a region in the North Central Pacific. 


B. DATA FORMATS AND PROCESSING 

This section describes the various formats in which the SSM/I and GEOSAT 
data are found. In addition, the process in which raw SSM/I information is decoded 
and then combined with the preprocessed GEOSAT wind speeds to form the final data 


base will be discussed. 


1. SSM/F Data Base 
There are four forms of data which are on file for each orbit made by the 
DSMP satellite. All are encoded in hexadecimal format and must be decoded to 
decimal before processing. 
The Orbit Summary File (SUM) contains information on the latitude and 
longitude of each end of the sensor swath for entire satellite orbit. Included with this 
is a history of sensor health, 1.e., whether or not any of the sensors on board 


experienced any unusual fluxes in radiation, including the hot calibration standard, 
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reference voltages, and the automatic gain control From this file, the exact time during 
which the DMSP orbit would coincide with a GEOSAT crossing could be determined. 
This fixed the amount of data needed to be downloaded to cover as much of the 
GEOSAT track as possible. 

Raw brightness temperatures for each sensor, are initially stored in a file known 
as the Sensor Data Record (SDR). An Antenna Patter Correction (APC) is applied 
to the brightness temperatures to correct for spurious energy received in the side lobes. 
Corrected brightness temperatures are then stored in the Temperature Data Record 
(TDR). It is in this format that brightness temperatures can be used to derive 
geophysical data such as precipitation rates, ice edge, ice age, water vapor, and wind 
speed. The geophysical data are then stored in the Environmental Data Record (EDR). 


This is the level of data used in this study. 


2. GEOSAT Data Base 

GEOSAT altimeter data is stored on board the satellite for approximately 
12 hours and then transmitted to a receiving station at the Johns Hopkins Applied 
Physics Laboratory (APL) where it is archived on tape. Three data products are 
produced from the raw data; The Sensor Data Record, the Naval Ocean Research and 
Development Activity (NORDA) Data Record, and the Waveform Data Record. The 
first two contain information about the measured altitude, significant wave height, wind 
speed, automatic gain control, and corrections for satellite and instrument errors. The 
last record contains waveform data used primarily for ice studies. 

FNOC stores GEOSAT data received from NORDA in half hour blocks 


which contain latitude, longitude, wind speed. and sea surface height information. 


oy: 


The data is pre-processed and does not contain the raw data necessary to verify the 
algorithms used to produce the wind speeds and corrections made to the data set. 
Reliance was placed on the ground truth comparison made with the GEOSAT wind 


speed calculations to establish the validity of the wind speed data. 


3. Processing GEOSAT Altimeter and SSM/I data sets 

Several Fortran programs were used to convert the hexadecimal format of 
the SSM/I data records to decimal. A program called COMPARE was written to 
combine the measurements from both the GEOSAT and SSM/I EDR data bases for 
each co-located satellite pass. SSM/I wind speeds measured within a 0.25 degree 
radius of a GEOSAT data point were averaged together and a mean value and standard 
deviation were obtained. A new data set was then produced which contained the 
GEOSAT wind speeds and the mean values of the SSM/I wind analysis from which 
correlation studies were made. A graphical representation of the two data sets were 
produced. The SSM/I EDR also supplied the rain, cloud water, and water vapor 
analysis necessary to determine the environmental conditions that existed during the 


wind speed measurements. 


C. DATA ANALYSIS AND RESULTS 

A comparison of all the compiled data sets for the SSM/I and GEOSAT altimeter 
wind speed analysis is seen in Fig. 6. The correlation for the compiled data sets is 
0.703. From this graph, four distinct regions are seen. The region "A" extends from 
2.5 to 9 m/s. Here, the SSM/I demonstrated good correlation with the GEOSAT 


altimeter. Wind speeds between 5 and 8 m/s represent compiled data from four 
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different orbits, however, wind speed from 2.5 to 5 m/s was derived from a single data 
set, which resulted in a region of sparse data. The SSM/I is limited to a minimum 
detectible wind speed of 3 m/s. This limitation is due to the lack of sensitivity to 
changes in emissivity for wind speeds between 0 and 3 m/s. The region greater than 
8 m/s demonstrated very poor correlation between the SSM/I and the GEOSAT 
altimeter. This area is divided in to three separate regimes. The region "B" 
corresponds to the bulk of data which is greater than 9 m/s. Here, the correlation 
appears to be good, but the SSM/I data are notably biased to the higher wind speeds. 
On either side of region "B", there are outlying points. Region "C" is identified by 
the extremely high values of wind speed of the SSM/I for GEOSAT winds from 8 to 
12 m/s. Region "D" is represented by the high values of the GEOSAT wind speeds 
relative to the SSM/I values from 11 to 14 m/s. 

As both systems have successfully verified wind speeds in the regime between 
8 to 15 m/s, the regions where the two measurements differed were examined in detail. 
The water vapor, cloud water content, and rain measured by SSM/I were used to 
determine if any interference due to these substances may have caused the bad 
correlation in the high wind speed areas. This was done for each individual data sets 
to determine if regional differences may have made a contribution. Data sets are 
identified by the corresponding SSM/I orbit number, latitude and longitude of the 
region, as well as the date and time difference between the SSM/I and GEOSAT co- 


located orbits. 
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Figure 6 SSMI versus GEOSAT wind speeds in meters per second. 
Line number 1 represents the actual correlation of 0.703 and line 2 
represents a reference correlation of 1.0. 


1. Orbit 9138, 27 March, 1989 


The scan area of the SSM/I encompassed a region bounded by the following 


coordinates: 


Latitude Longitude 
32.50°N 226.00°E 
32.50°N 242.20°E 
17.00°N 226.00°E 
17.00°N 242.20°E 
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The time difference between the GEOSAT and SSM/I passes is no greater than 0.61 
hours. The track of the GEOSAT through the SSM/I swath begins at 20°N,231.8°E 
and ends at 29.5°N,227.4°E. The wind speeds recorded in this region ranged from 
2.8 to 7.9 m/s. The correlation for this pass is 0.925 (Fig. 7). The difference in 
wind speed for the individual satellites demonstrates a maximum of 1.9 m/s (Fig. 8). 
The wind field, as observed by the SSM/I, is depicted in Fig. 9. The entire area of 
co-located orbits reveals a region of high water vapor content of approximately 60 
kg/m’ (Fig. 10). No cloud water or precipitation is present in this region. Since the 
correlation of the wind speed for orbit 9138 is so high, one could safely deduce that 
a water vapor content as high as 60 kg/m’ does not adversely affect the remote sensing 


of wind speeds by microwave instruments. 


2. Orbit 8934, 13 March, 1989 


The scan area of the SSM/I encompassed a region bounded by the following 


coordinates: 
Latitude Longitude 
57.94°N 2143325 
54.44°N 19259 a 
47.39°N Z1G97 E 
44.57°N 199.38°E 


The tume difference between the GEOSAT and SSM/I passes is within 1.01 hours. 


The track of co-located orbits begins at 39.5°N,217.8°E and ends at 27.8°N,211.5°E. 


36 


Soe ee GE OoAT WEN SPEED 


ORBIT 9138 


SOw/1 WINO SPEED (w/S) 





Q 2 4 6 8 10 
ALTIMETER WINO SPEED CWS) 


Figure 7 SSM/I vs altimeter wind speed. Line number | represents the 
actual correlation of 0.925 and line 2 represents a reference correlation 
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Figure 8 SSM/I and GEOSAT wind speed vs. latitude: orbit 9138. 





Figure 9 SSM/ wind speed field (m/s): orbit 9138. Line through field represents 
GEOSAT flight path. 
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Figure 10 SSM/I water vapor field (kg/m’); orbit 9138. The line through the field 
represents the GEOSAT flight path. A high correlation of wind speeds exist 
between the GEOSAT altimeter and SSM/I despite high water vapor concentration. 
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A break in the track was due to an interruption in communications during the transfer 
of data over a modem. Data were consequently lost between 31°N and 32.5°N 
latitudes. The statistics for this orbit revealed a correlation between the GEOSAT 
altimeter and SSM/I wind speeds of 0.616. This data set spans a range of wind speeds 
from 5 to 9.5 m/s (Fig. 11). The data points appear to be slightly biased toward the 
higher wind speeds for the GEOSAT altimeter. The plot of wind speed versus latitude 
(Fig. 12) reveals that both satellites follow the trend of increasing and decreasing wind 
speeds. The largest difference in wind speeds (approximately 2.5 m/s) appears at 
29.2°N, 32.5°N, 35.5°N, and 38°N latitudes. The images of wind speed (Fig. 13), 
water vapor (Fig. 14), and precipitation (Fig. 15) show that the track of the co-located 
orbit does not pass through any interfering atmospheric substances and that there were 
no unusual gradients in the vicinity. Overall, orbit 8934 demonstrated good correlation 


between the GEOSAT altimeter and the SSM/I. 


3. Orbit 8058(7), 10 January, 1989 


The scan area of the SSM/I encompassed a region bounded by the following 


coordinates: 
Latitude Longitude 
Bhetiein | Zivoe BS 
57.84°N 03 
51.88°N 223.64°E 
48.89°N 205-002 
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Figure 11 SSM/I vs. GEOSAT WIND SPEED: ORBIT 
8934. Line 1 represents the linear regression line. Line 
2 1S da reterence stove or 1.0. 
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Figure 12 SSM/I and GEOSAT winds vs. latitude: orbit 
8934. 
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Figure 13° SSM/I wind speed field (m/s): orbit 8934. Line through the wind 
field represents the GEOSAT flight path. 
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Figure 14 SSM/1 water vapor field (kg/m’); orbit 8934. Line through field 
represents the GEOSAT flight path. 
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Figure 15> SSM/I precipitation field (mm/hr): orbit 8934. Line through field 
represents the GEOSAT flight path. 





The time difference between the GEOSAT and SSM/I passes is no greater than 1.10 
hours. The track of the GEOSAT through the SSM/I swath begins at 49°N,219°E and 
altimeter for this particular orbit. By inspection, the correlation between the two 
instruments is poor. The calculated correlation is 0.208. Fig. 17 wlustrates the 
variation in wind speeds going from latitude 49°N to 57.5°N. Beginning at the most 
southem point, the difference between the SSM/I and GEOSAT altimeter grow from 
approximately 1 m/s to as much as 3.5 m/s. The corresponding SSM/I wind field 1s 
represented in fig. 18. Examination of the SSM/I water vapor field (Fig. 19), shows 
the orbit tracks cross through a gradient of water vapor which changes from 
approximately 10 to 30 kg/m’ over a distance of about 50 nm. The track continues 
through an area of saturated water vapor, i.e., greater than 60 kg/m’ and then exits into 
a region of approximately 30 kg/m’. The crossing of the gradient at 51°N coincides 
with the rapid increase in the SSM/I wind speed, while the altimeter remains relatively 
constant. Examination of the SSM/I precipitation field (Fig. 20) reveals precipitation 
rates from | to 4 mm/hr at 51°N and from 53.5°N to 56°N, along the co-located orbit 
track. From 56°N to the end of the track, cloud water is present up to 0.125 kg/m’ 
(Fig. 21). This suggests that in regions of cloud water and precipitation, SSM/I winds 


are greater than GEOSAT values. 
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Figure 16 SSM/I vs. GEOSAT wind speed: orbit 80587. Line |] 
represents an actual correlation of 0.208, while line 2 represents the 
reference correlation of 1.0. 
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Figure 17 SSM/I and GEOSAT wind speed vs. latitude: orbit 80587. 
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Figure 18 SSM/I wind field (m/s), orbit 8058(7). Line through field represents the 
GEOSAT flight path. 
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Figure 19 SSM/ water vapor field (kg/m’): orbit 8058(7). Line through field 
represents the GEOSAT flight path. 


48 


ee we £8: 

Sur ey: 
“2Pa'e 2h 0 PP 

s% Fak 0 eee wre . e 

PO ah al 

ee 

” Gira < 


Th on 
a} 


Fi . 
MMH: 


’ A pee 7 " FR 3, 
. Ky See htt * “Rea 2 } om nee Pn? 
BU Rigen Pc oat agi 
+ as eae eee, FH ¢é t"e a ° q's 
£ ‘ > eee s : 
RS key. gh F. a 
- ® ee 
é Phy “5 aod . we 
J a. © 
ated gt aan Ts ar = 
=e, 
Bogor bgh lyon 
Drea re ane oF 'E: 
PES ee 


ove 


, wit agree 2 ge % 
7 «a Ww 

on ae 7 
a i s 

s: &,' oe 

ye Be Sa) oe 5 


| ago ~ or 
fan. rh] 


: bs tis 4 wae | gees aig a te 2 
a! t - + ye sp Sage ol@s yA 
: « ea — 4.< r P-% 2 an 


S 





Line through wind 


field (mm/hr): arhit 8O5R8(7). 
49 


ipitation 


EFOSAT flight path. 


eure 20 SSM prec 
Id represents the G 


fie 


s 
‘t 





Tine throngh field 


arbit 8OS8(7). 


° 
€ 


SM/T cloud water field (kg/m’) 


represents the GEOSAT flight path. 


Figure 21S 


50 


4. Orbit 8058(6), 10 January, 1989 
The scan area of the SSM/I encompassed a region bounded by the following 


coordinates: 


Latitude Longitude 
51.88°N eo BLO4: E 
48.89°N 205.0 
42.14°N page) 0) ad 3 
39.60°N AAI 


The time difference between the GEOSAT and SSM/I passes is within 1.00 hours. 
The track of the GEOSAT through the SSM/I swath begins at 40.8°N,224.6°E and 
ends at 49.5°N,218.6°E. The corresponding graph of SSM/I vs GEOSAT wind fields 
is presented in Fig. 22. This data set spans a range of wind speeds from 7 to 12 m/s. 
The correlation in this graph is good below 8 m/s and then the SSM/I appears to read 
slightly higher then the GEOSAT altimeter. Overall, the correlation is calculated to 
be 0.935. A graph of wind speed versus latitude (Fig. 23) shows that both instruments 
were within 1 m/s from 41°N to about 46.5°N and afterwards the difference began to 
Open to about 2 m/s near 48.5°N. Examination of the SSM/I wind field (fig. 24) 
Shows that the location where the difference in wind speed begins to grow occurs 
around the same area that the wind speeds begin to increase rapidly. The SSM/I water 
vapor field (Fig. 25) does not demonstrate any significant gradient which would 


interfere with the measurement of wind speed or any area on increased water vapor. 


Sil 


SSM/ 1 VS ALTIMETER WIND SPEED 


ORBIT 80386 


SSw/1 WIND SPEED CWS) 





6 00 8 00 10 00 12 00 


GEOSAT WINO SPEEO (M/S) 


Figure 22 SSM/I vs. GEOSAT wind speed: orbit 80586. Line 1 
represents an actual correlation of 0.935, while line 2 represents a 
reference correlation of 1.0. 
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Figure 23 SSM/I vs. GEOSAT winds vs. latitude: orbit 8058(6). 
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Figure 28 SSMA/ water vapor field (kg/m’) 
represents the GEOSAT flight path. 
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There is, however, a gradual increase in water vapor content in the northeast quadrant 
of the field from 15 to 30 kg/m’. The precipitation field (Fig. 26) shows precipitation 
to the northeast of the GEOSAT-SSM/I track, but again, there is no evidence of 
precipitation occurring along the track. Overall, there appeared to be no interfering 


atmospheric substances 1n the region of the co-located orbits. 


5. Orbit 7925, 2 January, 1989 
The scan area of the SSM/I encompassed a region bounded by the following 


coordinates: 


Latitude Longitude 
79.98°N 33599-E 
70.99°N 296.68°E 
S2.c7 IN 348.24°E 
49.25°N 329 42°E 


The time difference between the GEOSAT and SSM/I passes is less than 0.50 hours. 
The track of the GEOSAT through the SSM/I swath begins at 54.6°N,343°E and 
terminates at 65.4°N,324.4°E. The plot of SSM/I versus GEOSAT wind speeds (Fig. 
27) demonstrates a tremendous amount of scatter in the data. Wund speeds range from 
7.5 m/s to 20.5 m/s. The correlation calculated for this data set 1s -0.301. Looking at 
the plot of wind speed versus latitude (Fig. 28), one can see that there are two distinct 
regions where the difference in wind speeds are greater than 3 m/s. This occurs 
between 57.5°N and 60°N and again from 63.5°N until 65.3°N. The SSM/I wind field 


(Fig. 29) exhibits a range of wind speeds from 5 to 31 m/s. The highest wind speeds 
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occur along a band extending from the south central region of the field toward the 
coast of iceland. High wind speeds are also observed along the coast of greenland. 
The SSM/I water vapor field (Fig. 30) exhibits a high concentration of water vapor in 
the southeast quadrant of the SSM/I and GEOSAT orbits. The large white area in the 
water vapor image represents values higher than 60 kg/m’. An examination of the 
precipitation field (Fig. 31) discloses a major region of precipitation occurring along 
the track as well. The location of the precipitation coincides with the same area of 
Fig. 28 that reveals a large difference in wind speed between the SSM/I and GEOSAT 
altrmeter (57.5°N to 60.0°N). The extent of precipitation appears to cover too small 
a region to be considered the only reason for the disparity between the GEOSAT and 
SSM/I wind speeds. The SSM/I cloud water field (Fig. 32) shows that this entire 
region is covered by high cloud water. From this evidence, the cloud water and 
precipitation together have apparently caused the large difference in wind speed values 
between the SSM/I and GEOSAT altimeter. 

A second area of concern lies between 63.5°N and 65.3°N (Fig. 28). This area 
is very close to the coast of Greenland, approximately 120 km at its nearest point. 
The altrmeter wind speed measurement steadily increases as it approaches Greenland, 
while the SSM/I wind speed measurements decrease. This phenomena may be due to 
the closeness of the land mass and side lobe effects on the SSM/I. Criteria for 
nearmess to land is not clear at this time. Chelton and McCabe [Ref. 18] and Wentz 


et al [Ref. 19] consider only those regions where microwave remote sensing is done 
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Figure 27. SSM/I vs. GEOSAT wind speed: orbit 7925. Line 1 
represents an actual correlation of -0.301, while line 2 represents a 


reference correlation of 1.0. 
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Figure 28 SSM/I and GEOSAT wind speeds vs. latitude. 
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Figure 29> SSM/T wind speed field (m/s): orbit 7925. Line represents the 
GEOSAT flight path. 
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Figure 30 SSM/T water vapor fteld (kg/m’): orbit 7925. Line through field 
represents the GEOSAT flight path. 
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from 200 to 800 km from the nearest Jand mass. They sight erratic behavior due to 
side lobes and the change in emissivity, when making microwave measurements of the 


Open ocean are made close to land. 


VIL. CONCLUSIONS 


To assess the relationship between the wind speeds as measured by the 
GEOSAT altimeter and the SSM/I, the correlation (R), R’, slope and intercept were 
calculated. Table one contains the results of the analysis of each data set. 

The value of R* represents the proportion of total variance in the values of the 
SSM/I wind speed data that can be accounted for by a linear relationship to the 
GEOSAT wind speeds. In orbits 9138 and 8058(6), more than 90% of the variance 
is explained by this linear relationship. Orbits 7925 and 8058(7) demonstrate almost 
no correlation between the SSM/I and altimeter wind speeds. 

The intercept gives the bias between wind speed measurements. Orbit 9138 
and 8058(6) exhibit the smallest bias, approximately 1 m/s. The accuracy goal of 

Table 1. Each data set is identified by its orbit number along with its 


associated correlation (R), R’, intercept, slope, and atmospheric parameters 
present during the wind speed measurement. 


ORBIT SLOPE INTERCEPT COMMENTS 


9138 io: 1.284 -1.297 High water vapor 
concentration 


8934 0.644 6.889 Clear Atmosphere 
8058(7) 10.599 Precipitation, cloud water, and 
water vapor 


8058(6) ACJ) -1.800 Clear Atmosphere 


7925 16.980 Precipitation, cloud water, and 
water vapor 





the GEOSAT altimeter and the SSM/I has been 2 m/s. Both orbits 9138 and 
8058(6) are within that goal. Orbits 7925, 8058(7), and 8934 have a bias in wind 
speed greater than 2 m/s. The slope is also a good indicator of the relationship 
between the SSM/I and GEOSAT altimeter wind speeds. A perfect relationship 
between wind speed measurements would give a slope of 1.0. Orbits 9138 and 
8058(6) have a slope very close to one, while 7925 and 8058(7) have a slope near 
0.0. 

There are four conclusions that are made from the analysed data in this paper: 

(1) The data sets which did not contain cloud water or precipitation along the 
co-located SSM/I and GEOSAT orbits (orbits 9138, 8934, and 8058(6)), 
demonstrated a good correlation between the two systems. 

(2) Although a large amount of water vapor was present along the entire track 
of orbit 9138, there appeared to be no evidence of bias with the measured wind 
speeds from either instrument. Since the SSM/I corrects for the presence of water 
vapor in the atmosphere, it is obvious, from the high correlation between the 
GEOSAT altimeter and the SSM/I, that the SSM/I water vapor correction works 
well. The altimeter was not corrected for water vapor, suggesting that water vapor 
corrections on the altumeter were small. 

(3) The data sets which contained cloud water and precipitation (orbits 8058(7) 
and 7925), displayed large differences in recorded wind speeds between the SSM/I 
and GEOSAT altimeter. In these areas, the SSM/I wind speed would increase, 
while the altimeter wind speed decreased relative to one another and, consequently, 


very poor correlation resulted. There are several explanations for the difference in 
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wind speed relative to each other experienced by both sensors. The increased 
emissivity associated with precipitation, or cloud water, is the most likely cause of 
increased wind speeds experienced by the SSM/I, by introducing a "warm " 
atmospheric term. The GEOSAT may experience a decrease in wind speed at the 
same location. This may have been a result of a decrease in ocean surface 
roughness. The reduction in surface roughness is caused by rain drops impacting on 
the ocean surface and consequently attenuating the surface gravity waves [Ref. 21] . 
The reduction in surface roughness would result in a larger backscatter cross section 
for the GEOSAT altimeter. This reduced backscatter is interpreted as a decrease in 
wind speed, which is what is observed in the area of precipitation in orbits 7925 
and 8058(7). An alternative explanation for the low GEOSAT wind speed in the 
area of precipitation is that the measurement was correct and that the SSM/I had a 
bias to high wind speeds in the presence of rain. 

(4) Erratic behavior is seen in microwave measurements of the open ocean 
near a land mass. This is evident in orbit 7925 when measurements of wind speed 
and atmospheric properties were made within 120 km of the Greenland coastline. 
The resulting data from the coastal areas was highly erratic. Based on only one 
data set near a coastal feature, it appears that wind speed measured within 200 km 
of a land mass is not reliable. This is most likely due to the radiance detected by 
antenna sidelobes which mask the actual measurement of radiation from the main 
lobe of the antenna. 


The wind speed verification studies of the GEOSAT altimeter and the SSM/I 


demonstrate that both active and passive microwave systems measure similar values 
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of wind speed over the ocean. These instruments are comparable in open ocean 
regions at least 200 km from land and free of cloud water and precipitation. 
Utilization of GEOSAT altimeter and SSM/I wind speed data in a numerical 
forecasting model requires that a system for screening wind data for cloud water, 


precipitation, and proximity to land masses be used. 
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